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^ ' Abstract. We study an incommensurate long-range order induced by an 

Q ' external magnetic field in a quasi-one-dimensional bond-alternating spin system, 

^ ' F5PNN, focusing on the role of the frustrating interaction which can be enhanced 

by a high-pressure effect. On the basis of the density matrix renormalization 

group analysis of a microscopic model for F5PNN, we present several H-T 

phase diagrams for typical parameters of the frustrating next-nearest-neighbour 

^ ' coupling and the interchain interaction, and then discuss how the field-induced 

^^ ' incommensurate order develops by the frustration effect in such phase diagrams. 

^y\ A magnetization plateau at half the saturation moment is also mentioned. 
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1. Introduction 

Quasi-one-dimensional quantum spin systems exhibit several interesting field-induced 
phenomena. In a gapped spin chain an external magnetic field exhibits a quantum 
phase transition to the gapless Tomonaga-Luttinger liquid phase at some critical 
field Hci P EJ . The Tomonaga-Luttinger liquid phase basically continues up to the 
saturation field Hc2 and is characterized by the power-law decay of the spin correlation 
function: 

{SSS:)r^i-iyr-^^, (1) 

(S^S^) - m^ - cos{2kFr)r-'^' , (2) 

where m is the magnetization along the field H \\ z, and the Fermi wave number 
kp depends on the magnetization. Since these exponents satisfy rj^ < rjz in usual 
antiferromagnets, the commensurate antiferromagnetic long-range order perpendicular 
to H should appear in the presence of interchain interactions. In fact, such 
a field-induced antiferromagnetic order has been theoretically proposed [HI and 
experimentally observed in some quasi-one-dimensional gapped systems |2|2]. Here 
it should be noted that the ordered phase can be also interpreted as a result of the 
magnon Bose- Einstein condensation |5|. 

In some one-dimensional (ID) frustrated systems, however, the 77-inversion, 
namely rjr^ > rj^, has been predicted for some intermediate strength of the magnetic 
field between Hd and Hc2 [3 El 111 ^1 Qj- Indeed, such ?;- inversion has been 
suggested experimentally by the NMR relaxation rate 1/Ti of an organic S = 1/2 
bond-alternating chain, pentafulorophenyl nitronyl nitroxide (F5PNN), in a certain 
temperature region where a power-law behaviour of l/Ti can be observed |12l I18j . 
A more interesting point associated with the ry-inversion is that the field-induced 
long-range incommensurate order can appear, refiecting the enhancement of the 
incommensurate correlation ^31- Although the recent specific heat measurement |15l 
I16| indicates that the field-induced commensurate antiferromagnetic order may be 
realized for F5PNN, the NMR experiments are still suggesting that an enhancement 
of the incommensurate correlation in F5PNN. Thus, it seems that F5PNN is located 
in the vicinity of the boundary between the commensurate order phase and the 
incommensurate order phase. In this paper, we discuss the possibility of realizing 
the field-induced incommensurate order in F5PNN, putting a special emphasize on 
a high-pressure effect which can enhance the frustrating interaction. We examine 
how the incommensurate order appears in the phase diagram on the basis of density 
matrix renormalization group calculations of a microscopic model for F5PNN. We also 
mention a possibility for a magnetization plateau and coexistence of the commensurate 
and incommensurate orders. 

2. Model 

The experimentally observed ?7-inversion and the field-dependent crossover of the 
bond-alternating ratio J7] have indicated the presence of frustration in F5PNN. 
Thus, we consider the bond-alternating spin chain with the next-nearest-neighbour 
interaction in a magnetic field 

H = Ji 2_^ S2J ■ S2J+I + J2 2_^ ^^j + l ■ '^2j-|-2 (3) 
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+ J'J2^2rS2j+2-HY,S'-, (4) 

J 3 

as a realistic model for F5PNN, where Ji and J2 {^ Ji) denote the alternating nearest 
neighbour interactions, and J' means the next-nearest neighbour coupling. In what 
follows, Ji is set to unity for simplicity. For Hd < H < Hc2, the model is basically in 
the Tomonaga-Luttinger liquid phase with the Fermi wave number kp = (1/2 — m)?:. 

3. Mechanism of the field-induced incomniensurate order 

Let us briefly explain the mechanism of the field-induced incomniensurate order. If 
J' is sufficiently large, the ?7-inversion occurs around m ~ TOs/2, where m^ is the 
saturated magnetization |S1 |H1 ^3 . At the magnetic field in the yy-inversion region, the 
dominant spin correlation should not be the usual commensurate antiferromagnetic 
one perpendicular to 7J, but the incommensurate one along the _ff-direction. In the 
presence of interchain interactions, the long-range order corresponding to the dominant 
spin correlation is expected, implying that a long-range incommensurate order in the 
H direction can be realized for the region rj^ > rjz- 

In the present study, the dominant spin correlation is determined by two 
staggered susceptibilities X-L and X|| corresponding to the spin correlations Q) and 
(0 , respectively. Low temperature behaviours of the susceptibilities are characterized 
by the exponents 7]x and r]z |1U) : 

Xx(T) ^ C^T-^^-"-) and x\\{T) ~ G\\T-'^^-'''\ (5) 

where Cj^ and Cii arc nonunivcrsal coefficients. Equation |(SJ) shows that the more 
dominant spin correlation corresponds to the more strongly divergent susceptibility. 
When the 77-inversion occurs, for example, X|| shows stronger divergence than x±- 

For quantitative treatment of the field-induced order and of the accompanied 
phase transition, we use the above susceptibilities. Based on the mean-field 
approximation for the interchain interaction Jint |18| . the critical temperature Tc is 
given by 

X^{T,) = (zJint)-\ (6) 

where 7 is _L or ||, and z the number of adjacent chains. 

On the basis of equations © and © , we can classify possible types of the field- 
induced orders; we consider the following cases, which are depicted in figure 1. The 
first one is that -q^ < rjz [figure Q] (a)] ^U]. In this case, since X-L strongly diverges, 
X-L > X\\ is always realized. As T decreases, X-L first satisfies the relation © and 
the commensurate order appears below T^. By contrast, if the ry- inversion occurs 
(i.e. rjx > r]z), x± < X\\ is naively expected. However, we should pay attention to 
the magnitude of the non-universal coefficients C_l and C||. We have found that, 
even for rj^ > rjz, C_l is larger than C||, and thus xJ- can be larger than X|| at high 
temperatures, although x\\ finally exceeds X-L at T = T* . If T^ > T*, therefore, x± 
satisfies equation Q at a higher temperature and the commensurate order appears 
[see figure n(b)]. In such a case, another phase transition is expected to occur from 
the commensurate phase to the incommensurate phase at T = T*. The previous 
DMRG analysis of the free energy has illustrated that the transition actually occurs 
when zJint is not so large ^3. For Tc < T* [figure ^(c)], X|| satisfies equation 
first and only the incommensurate order appears . 
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Figure 1. Schematic temperature dependence of x± a-nd Xll ^^ three cases. At 
T = Tc, the larger susceptibility satisfies the relation J^) . At T = T* in (b) and 
(c), x± is equal to xil- 



4. Phase diagrams 

By a numerical calculation of xi- ^-i^d x\\ ^^r a single chain and the mean-field 
treatment described above, we have obtained H-T phase diagrams. For the calculation 
at finite temperature in the thermodynamic limit, we use the finite temperature 
DMRG |21[ 121] . The DMRG calculation has been carried out with the maximum 
number of the retained bases m = 64 and the truncation error is at most the order of 
10~^. Technical details to obtain x± and x\\ are described in reference .14 . 

In figure El we show the phase diagrams for typical three parameters of the 
frustrating coupling: (a) J' — 0.05 (small J' case), (b) J' — 0.1 (intermediate J' case), 
and (c) J' — 0.15 (large J' case). The interchain interaction is fixed to zJint = 1/12, 
and J2 is fixed to a realistic value 0.45 for F5PNN ^J. In the figure, C and IC 
denote the commensurate antiferromagnetic-order phase and the incommensurate one, 
respectively. 

4.I. case (a) 

The phase diagram in the case (a) is typical for spin-gapped quasi-lD spin systems, 
where it should be remarked that the shape of the phase boundary between the C 
phase and the disordered phase looks like a semicircle. Since x± is always larger than 
X\\ as shown in figuren(a), the IC phase never appears. 



4-2. case(h) 

In the case (b), the phase diagram is clearly different from the case (a) in the following 
two points. The first one is the appearance of the IC phase in a very low temperature 
region. At H = 1.23, where the magnetization m at T = is equal to TOs/2, the 
ry-inversion occurs and thus xii shows slightly stronger divergence than that of x±- 
Thus the IC phase can appear around H = 1.23 in the quite low temperature region 
(T < T*). However, in the high temperature region (T > T*), xa. is larger than xy, 
resulting in the transition between the commensurate order phase and the disordered 
phase at T = Tc. These behaviours of the staggered susceptibilities correspond to 
figure n(b). Another important point is that the phase boundary between the C 
phase and the disordered phase is deformed from the semicircle-shape particularly 
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Figure 2. H — T phase diagram for various J' (J' = 0.05, 0.1 and 0.15). C 
and IC denote the commensurate antiferromagnetic and incommensurate order 
phases, respectively. Filled circles denote the transition temperature between the 
IC phase and other phases, and open circles indicate that between the C phase 
and the disordered phase. The interchain interaction is fixed to zJi^t = 1/12. 



around H ~ 1.23. This is because the transverse fluctuation is suppressed by the 
frustration effect even in T > T*, resulting in the lowering of Tc- 

4-3. case (c) 

In the case (c), the large J' enhances the incommensurate correlation and extends the 
region of the IC phase. Around H — 1.27, xj_ behaves as in figure^] (c) and thus we 
can observe the direct transition from the disordered phase to the IC phase, which 
divides the C phase into the \avj-H part and the high-iJ part. 



5. Discussion 

In the previous work TT, the authors estimated J' of F5PNN as ~ 0.05(= 0.3K), based 
on a qualitative analysis of the magnetization curve. This value is slightly smaller than 
the critical value (~ 0.08) for the ry-inversion or the emergence of the incommensurate 
order (see figure 13) . We however think that high-pressure experiments can realize 
the incommensurate phase within realistic experimental situations. For example, a 
hydrostatic pressure experiment reveals a drastic change of the bond-alternation ratio 
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Figure 3. J2-J' phase diagram of the dominant correlations at rra = ms/2 
at zero temperature, which is equivalent to the 1/2 plateau phase diagram in 
reference 1261 . The filled circle corresponds to the parameters of F5PNN. 



in powder samples of F5PNN ,23,. In addition, a recent experiment clarifies that 
specific heat of the power samples has slightly different temperature dependence from 
that of a single crystal of F5PNN, which can be attributed to an effective pressure 
caused by an experimentally technical reason .24 . These experimental results strongly 
indicate that the properties of F5PNN are quite sensitive to the pressure. Thus it is 
worthwhile to discuss precisely how the incommensurate order phase develops in the 
H-T phase diagram near the phase boundary of the 77-inversion. 

As in figure O we can readily consider possible two paths such that the couplings 
of F5PNN exceeds the critical value for the incommensurate order: increasing the next- 
nearest neighbour interaction J' (arrow a) and/or decreasing the alternating coupling 
J2 (arrow /3). Here, we chiefly discuss the J'-dependence of the H-T phase diagrams 
along the vertical arrow direction in figure 13 assuming that Jint is independent of the 
pressure. We can then expect that the features of H-T phase diagrams should change 
as figure 121 (a) -^ figure [3 (b) -^ figure |21 (c). If J' is sufficiently enhanced by the 
pressure, the direct phase transition from the disordered phase to the IC phase would 
be observed in a middle range of the magnetic field as in the case (c). As mentioned 
in the previous section, however, a notable point is that, even in the case (b), we may 
find the deformation of the phase boundary between the disordered phase and the C 
phase. From the experimental view, such a deformation of the phase boundary can 
be an important signal of a development of the IC order induced by the frustration, 
though it may be difficult to capture the IC phase at a very low temperature. 

Next, we discuss how the changes of the exchange interactions can reflect on a 
measurement of the specific heat at the zero magnetic field, which is quite helpful in 
analysis of the experimental results, because it is not an easy work to obtain a complete 
H-T phase diagram in actual experiments [151 I16| . Figure 0] shows temperature 
dependence of the specific heat for various values of J' and J2 along the arrows in 
figure 13 When the frustrating interaction J' increases along the a direction, it is 
found that the pronounced peak around T = 0.3 sharpens. We note that this can be 
attributed to the increase of the density of state caused by the narrowing of the S = I 
magnon band width |25| . On the other hand, such a band narrowing also occurs, as J2 
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Figure 4. Temperature dependence of the specific iieat for various values of 
excliange parameters. Each figure corresponds to the case where the exchange 
interactions are varied along the arrows a and /3 in figure ISl 




Figure 5. H-T phase diagram in case of large J' and large Jint (</' = 0.15 and 

zJint = 1/6). 



is reduced along the /3 direction. We can then see that the peak sharpens as well. We 
therefore think that a careful analysis of the shape of the specific heat peak at i:f = 
provides an important information for the enhancement of the IC correlation in the 
intermediate magnetic field region, although the pressure effect subtly influences the 
both of J' and J2 in the actual situation. 

Here, we should make a comment on the pressure effect for the interchain 
interaction. If Jint is also enhanced by the pressure, the direct phase transition from 
the disordered phase to the IC phase may be suppressed within the interchain mean- 
field analysis as shown in figure |31 where the interchain interaction is assumed to 
be unfrustrating. Nevertheless, the precise analysis of the crystal structure ^Tj of 
F5PNN has indicated that the interchain couplings are also frustrating, which would 
effectively decrease value of J-mt in the mean-field approximation. In order to reveal 
the role of the frustrating interchain interaction, of course, a precise analysis beyond 
the mean-field level is clearly required. 
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Figure 6. Magnetization curves for the realistic parameters J2 = 0.45 and 
J' = 0.15 at several temperatures. 



Finally, we want to discuss a magnetization plateau at the half of the full- moment. 
As in figure 13 the 1/2 plateau is closely related to the mechanism of the 77-inversion; 
the criterion for the ?;-inversion is equivalent to the one for the plateau formation jl4| . 
The level spectroscopy analysis [201 on the model l@J has indicated the 1/2 plateau 
should appear at half the saturated magnetization for Ji = 0.45 and J' = 0.15. 
However, we think that it may be difficult to observed the clear evidence of the 1/2 
plateau of F5PNN even with a high pressure experiment. This is because the width of 
the plateau is still very narrow for these couplings, and thus the thermal fluctuation 
smears out the plateau even at temperatures where we can observe the magnetization. 
In figure El we show the magnetization curves for J2 = 0.45 and J' = 0.15 obtained 
by the finite temperature DMRG, illustrating that the plateau is easily smeared out 
by the thermal effect. We think that, as for F5PNN, it is rather hopeful to verify the 
presence of the IC order than to make a direct observation of the 1/2 plateau. 

To summarize, we have discussed field-induced IC order for the bond-alternating 
spin chain with the next-nearest interaction, which is relevant to the high pressure 
experiments of F5PNN. In particular, we have analyzed how the various physical 
quantities behaves near the phase boundary, and then concluded that the precise 
analysis for the phase boundary and the zero-field specific heat can detect a signal 
of the field-induced IC order. Further experimental studies on F5PNN are quite 
interesting. In addition, it is also a challenging issue to search another spin chains 
exhibiting the field-induced incommensurate order. A recently observed new phase of 
NDMAP in the magnetic field might be a possible candidate j27| . 

The field-induced commensurate order in a quantum spin system is often 
described as the Bose Einstein condensation (BEC) of a magnon gas. The present 
IC order may be viewed as a result of a "charge order" of the magnon. As mentioned 
in reference j2()| the transition between the field-induced incommensurate order and 
the usual commensurate antiferromagnetic one is of the first order. Then we can expect 
the coexistence of the two orders at the phase boundary. Recently, such a coexistence 
of the two different orders is actually observed in solid ^He, which is called "super- 
solid" state I2H1. A connection to such an exotic state is also an interesting problem 
both from experimental and theoretical viewpoints. 
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